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ABSTRACT 

Early phase optical spectra of aspherical jet-like supernovae (SNe) are presented. We focus on energetic 
core-collapse SNe, or hypernovae. Based on hydrodynamic and nucleo synthetic models, radiative transfer in 
SN atmosphere is solved with a multi-dimensional Monte-Carlo radiative transfer code, SAMURAI. Since the 
luminosity is boosted in the jet direction, the temperature there is higher than in the equatorial plane by ~ 
2,000 K. This causes anisotropic ionization in the ejecta. Emergent spectra are different depending on viewing 
angle, reflecting both aspherical abundance distribution and anisotropic ionization. Spectra computed with an 
aspherical explosion model with kinetic energy 20 x 10 51 ergs are compatible with those of the Type Ic SN 
1998bw if ~ 10-20% of the synthesized metals are mixed out to higher velocities. The simulations enable us 
to predict the properties of off-axis hypernovae. Even if an aspherical hypernova explosion is observed from 
the side, it should show hypernova-like spectra but with some differences in the line velocity, the width of the 
Fe absorptions and the strength of the Na I line. 

Subject headings: supernovae: general — supernovae: individual (SN 1998bw) — radiative transfer 
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1. INTRODUCTION 

The connection between the long gamma-ray bursts 
(GRBs) and a special class of Type Ic supernovae (GRB-SNe) 
is now well established (Galama et al. 1998; Hjorth et al. 
2003; Stanek et al. 2003; Malesani et al. 2004; Pian et al. 
2006). Since GRBs are induced by relativistic jets, GRB-SNe 
are also expected to be aspherical. There is also increasing 
evidence that core-collapse SNe are not spherically symmet- 
ric, coming from, e.g., the detection of polarization in several 
SNe (e.g., Wang et al. 2001; Kawabata et al. 2002; Leonard 
et al. 2006) and late-time spectroscopy (Mazzali et al. 2005). 

However, the progenitors, the explosion mechanisms and 
the origins of diverse properties of GRB-SNe are still not un- 
derstood (e.g., Nomoto et al. 2007). In order to answer such 
questions, it is crucial that the properties of the explosions 
(e.g., the mass and kinetic energy of the ejecta and the asym- 
metry of the explosion) are accurately derived from observa- 
tions. For this purpose, we should know the observational 
properties of multi-dimensional SN explosions, and to what 
degree they are affected by orientation effects. 

This is an area that has not been studied in great depth, de- 
spite its importance. Although some works have addressed 
the properties of the light curve (LC) or the late phase spec- 
tra of aspherical SNe (e.g., Hoflich et al. 1999; Maeda et al. 
2002; Maeda et al. 2006a, c; Sim 2006), only a few studies 
of early phase spectra (t < 50 days, where t is the time since 
the explosion) have been performed (e.g., Hoflich et al. 1996 
for core-collapse SNe, Thomas et al. 2002; Kasen et al. 2004; 
Tanaka et al. 2006 for Type la SNe). 

SNe that associate with GRBs are thought to be highly en- 
ergetic explosions, i.e., hypernovae (here defined as SNe with 
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ejecta kinetic energy E$i = £ K /10 51 ergs > 10; e.g., Nomoto 
et al. 2006) indicated by broad line features in their early 
phase spectra. However, this suggestion was based on anal- 
yses that assumed spherical symmetry. No realistic multi- 
dimensional explosion models have been verified against the 
observed early phase spectra. 

We report the first study of early phase spectra of aspherical 
hypernova models using a multi-dimensional radiative trans- 
fer code. The synthetic spectra are compared with observed 
spectra of SN 1998bw, and implications for off-axis hyper- 
novae are discussed. 

2. MODELS 

We use the results of multi-dimensional hydrodynamic and 
nucleosynthetic calculations for SN 1998bw (Maeda et al. 
2002) as input density and element distributions. In the 
present simulations, 16 elements are included, i.e., H, He, C, 
N, O, Na, Mg, Si, S, Ti, Cr, Ca, Ti, Fe, Co and Ni. Since 
the original models used a He star as a progenitor, we simply 
replace the abundance of the He layer with that of the C+O 
layer. In the hydrodynamic model, energy is deposited as- 
pherically, with more energy in the jet direction (z-axis). As 
a result, 56 Ni is preferentially synthesized along this direction 
(see Fig. 3 of Maeda et al. 2002). In this Letter, an aspheri- 
cal model with E$i = 20 (A20) and two spherical models with 
£51 = 20 and 50 (F20 and F50, respectively) are studied (Table 
1). They are constructed based on the models with E51 = 10 
(Maeda et al. 2002; Maeda et al. 2006c). 

Bolometric light curves (LCs) computed in 3D space by 
Maeda et al. (2006c) are also used as input for the spec- 
tral calculations. A common problem in modeling hypernova 
LCs is that a spherical model reproducing the early rise of the 
LC (£51 = 50 for SN 1998bw) declines more rapidly than the 
observed LC at t ~ 100- 150 days (Nakamura et al. 2001; 
Maeda et al. 2003). This problem can be solved by aspherical 
models with a polar view. In aspherical models, even with a 
lower kinetic energy (£51 = 10-20) than in the spherical case 
(£"51 = 50), which allows sufficient trapping of 7-rays at late 
times, the rapid rise of the LC can be reproduced because of 
the extended 56 Ni distribution (Maeda et al. 2006c). 



2 



Tanaka et al. 



TABLE l 

Summary of Explosion Models 



model name 




Esi 


M( 56 Ni) b 


P c 


A20 


10 


20 


0.39 


0.0 


F20 


10 


20 


0.31 


0.0 


F50 


10 


50 


0.40 


0.0 


A20p0.2 


10 


20 


0.39 


0.2 



a The mass of the ejecta (Mq) 
b The ejected 56 Ni mass (Af ) 
c The mixing fraction (see §5) 



3. METHOD 

In order to study the detailed properties of the radiation 
from aspherical SNe, we have unified a SupernovA Multi- 
dimensional RAdlative transfer code SAMURAI. SAMURAI 
is a combination of 3D codes adopting Monte-Carlo methods 
to compute the bolometric LC (Maeda et al. 2006c), and the 
spectra of SNe from early (Tanaka et al. 2006) to late phases 
(Maeda et al. 2006a). The LCs are computed as pioneered 
by Cappellaro et al. (1997). Our 3D LC code adopts individ- 
ual packet method as described by Lucy (2005) and addition- 
ally, includes multi-energy transfer for high energy photons 
(Maeda 2006b). 

The early phase spectra are calculated as snapshots in the 
optically thin atmosphere, using the results of the LC simu- 
lation as initial conditions. A sharply defined photosphere is 
assumed as an inner boundary for simplicity. The position of 
the inner boundary in each direction is determined by aver- 
aging the positions of the last scattering photon packets (see 
Fig. 3 of Maeda et al. 2006c). Thus, the inner boundary cor- 
responds to an optical depth of roughly unity. The luminosity 
at the inner boundary (Lm(9)) is also taken from the LC com- 
putation by integrating the energy of photon packets escaping 
from each bin. 

For the computation of ionization and excitation state in the 
atmosphere, the local physical process same as in the previ- 
ous ID code (Mazzali & Lucy 1993; Mazzali 2000), which is 
used to model spectra of SNe of various types (e.g., Mazzali 
et al. 1993), is adopted. Line scattering under the Sobolev 
approximation and electron scattering are taken into account. 
For line scattering, the effect of photon branching is included 
as in Lucy (1999). Beginning with a trial temperature struc- 
ture, a number of photon packets (typically ~ 10 8 - 10 9 ) are 
followed in 3D space, giving the flux in each mesh. Then the 
temperature structure is updated. After the temperature struc- 
ture converges, the spectra are obtained by counting the en- 
ergy (or wavelength) distribution of escaping photon packets 
in every direction. 

Since the background model is axially symmetric, it is not 
necessary to follow the propagation in the azimuthal direction 
explicitly. We use a 2D grid (r, 6), but still perform the trans- 
port of the photon packets in 3D, reproducing photon prop- 
agation in the azimuthal direction by projecting the photon 
back on the 2D plane of the calculation. The spatial grid has 
a typical mesh number (n n no) = (20, 80). Velocity is used as 
a proxy for the radial coordinate thanks to the homologous 
expansion of SN ejecta (r = vt). Details of the code will be 
presented elsewhere. 

4. RESULTS 

At t = 20 days, the photosphere of model A20 is mildly as- 
pherical with an axis ratio of 1 .4 between the polar and equa- 
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Fig. 1 . — Upper left: Computed temperature structure in the atmosphere of 
model A20 at 20 days after the explosion. Upper right, lower left: Ionization 
fractions of Fe II and Fe III. Lower right: Sobolev optical depth of the Fe II 
A5018. The vertical axis (z-axis) is the jet direction. 



torial directions. The emergent luminosity is also anisotropic, 
being brighter by a factor of 1.2 in the polar direction (Maeda 
et al. 2006c). In contrast, the local luminosity at the photo- 
sphere, L m (0), is highly aspherical, with an axis ratio of 6.3, 
tracing the aspherical 56 Ni distribution. At this epoch, the 
photosphere is still outside the region where explosive nucle- 
osynthesis occurs. The metal lines in the spectrum result from 
pre-explosion abundance in the C+O layer. 

The upper left panel of Figure 1 shows the calculated tem- 
perature structure. The boosted luminosity makes the tem- 
perature near the z-axis higher than in the equatorial plane 
by ~ 2000 K. This anisotropy of the temperature causes 
anisotropies of the ionization structure. The upper right and 
lower left panels of Figure 1 show the ionization fractions of 
Fe II and Fe III, respectively. Because of the high tempera- 
ture, Fe II is suppressed by a factor of more than 30 near the 
z-axis. As a result, the optical depth of the Fe II lines is larger 
in the equatorial plane (Fig. 1, lower right panel). 

The upper panel of Figure 2 shows the synthetic spectra at 
t = 20 days for models A20 (red, green and blue for = 0° (po- 
lar), 45° and 90° (equatorial), respectively), F20 (black) and 
F50 (gray). Here 6 is measured from the z-axis. All absorp- 
tion lines except for Si II A6355 are stronger for larger 9, i.e., 
for an equatorial view. This is because all species that have 
strong lines, i.e., O I, Si II, Ca II, Ti II, Cr II and Fe II, domi- 
nate near the equator but not near the z-axis as shown in Fig- 
ure 1 for the case of Fe. 5 . This also leads to a more effective 
flux-blocking in the near-UV in the side-viewed spectrum. 

The lower panel of Figure 2 shows the synthetic spectra 
at t = 30 days. At this epoch, the local luminosity and the 
temperature structure are still anisotropic, although the pho- 
tosphere is almost spherical (upper left panel of Fig. 3). Con- 
sequently, the distribution of ionization fractions is also as- 
pherical (see Fig. 3, upper right and lower left panels, for the 
case of Ca). However, the emergent spectra are not signifi- 
cantly different for different viewing angles (Fig. 2). This is 
the effect of the aspherical abundance distribution. The pho- 
tosphere at this epoch is located inside the region where heavy 

5 Si II A6355 is not strong at 6 = 90° because the temperature in the equa- 
torial plane is too low to activate the transition. 
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Fig. 2. — Upper: Synthetic spectra of model A20, F20 (black) and F50 
(gray) at 20 days after the explosion. The red, green and blue lines show the 
spectra of model A20 seen from 9 = 0° , 45° and 90° , respectively. Here is 
measured from the z-axis (jet direction). Lower: Same as upper panel but at 
30 days after the explosion. 



elements are synthesized in the explosion. Since nucleosyn- 
thesis occurs entirely near the polar direction in our model, the 
suppression of important ions near the z-axis is compensated 
by the larger abundance of the heavy elements (see Fig. 3, 
lower right panel) 6 . Only the Na I A5890 line is stronger for 
larger 6 because of the combined effect of the lower temper- 
ature that favors Na I and the higher abundance of Na, which 
is predominantly synthesized before the explosion. 

Given the strong anisotropy in L m (Q), the side- viewed spec- 
trum of model A20 consists mainly of photons escaping from 
the ejecta of 6 ~ 60°. Since physical quantities there (e.g., 
isotropic mass and kinetic energy of the ejecta in that direc- 
tion, and L m {6)) are similar to those of model F20, the spec- 
trum of Model F20 is similar to the polar- viewed spectrum of 
model A20. Model F50 shows a very strong Na I A5890 line 
owing to the low temperature. 

5. COMPARISON WITH SN 1998bw 

The spectra computed for the aspherical model A20 are 
compared with those of SN 1998bw in Figure 4. Because 
of the detection of GRB 980425, it is tempting to compare the 
polar- viewed spectrum (red lines) with that of SN 1998bw. 
Since the computed LC is brighter than the observation by a 
factor of 1.25 at t = 20 days (Maeda et al. 2006c), the syn- 
thetic spectra ait = 20 days are scaled after the computation 



6 The Fe II absorptions are broader at t = 30 days than at t 
because of more efficient absorption both near the photosphere (v 



20 days 
< 12,000 

km s _1 , a lower velocity than at t = 20 days) and at high velocity (v ~ 25, 000 
km s _1 , where the temperature is lower and the Fe II fraction is higher than at 
t = 20 days). 
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Fig. 3. — Upper left: Computed temperature structure in the atmosphere of 
model A20 at 30 days after the explosion. Upper right, lower left: Ionization 
fractions of Ca II and Ca III. Lower right: Sobolev optical depth of Ca II 
A8542. 



by this factor in Figure 4. 

At t = 20 days (Fig. 4, upper panel), the scaled synthetic 
spectrum of model A20 with = 0° does not show the broad 
absorption trough between 4000 and 5000A that is typical of 
hypernovae such as SN 1998bw. This is because in the model 
Fe-line absorption is not effective at high velocities, making 
it possible for photons to escape effectively near 4500A. Sim- 
ilarly, the O I-Ca II absorption at 7000 - 8000 A is very weak 
in Model A20. The Si line velocity in the model is, however, 
comparable to that in SN 1998bw. 

At t = 30 days (Fig. 4, lower panel), the Ca II and Fe II 
lines are somewhat stronger, although the O I-Ca II absorp- 
tion at 7000 - 8000 A is still narrower than in the observed 
spectrum. The peaks around 4000 and 4500 A are partially 
suppressed for a polar view by the high velocity absorption 
by the extended Fe near the jet, while they are strong in the 
synthetic spectrum for the equatorial view. The suppression 
of the peaks is similarly seen in the spectrum of SN 1998bw. 
The Na I A5890 line is very weak in SN 1998bw, in analogy 
with the spectrum with = 0°. 

The strengths of the Ca II and Fe II lines aU = 20 days can 
be increased if heavy elements synthesized in the explosion 
are mixed to outer layers. In SN explosions, Rayleigh - Tay- 
lor (R-T) instabilities are expected to occur 7 , which could 
deliver the newly synthesized elements to higher velocities. 
We simulate this possibility by introducing a parameter, the 
mixing fraction p. If the mass ejected inside a conical sec- 
tion of the ejecta is M e j(0), and Mi(0) is the mass of a certain 
element contained in the conical section, we simulate mixing 
by taking a fraction (p) of the mass of a newly synthesized 
element and distributing this mass homogeneously in the con- 
ical section. Thus, in the outer layers, the newly synthesized 
element will have an abundance X;(r, 0) = pMi(0)/M^{6) (+ 
pre-SN abundance). In the deeper region, the abundances are 
redistributed by a mixing down of the unburned material. 

In Figure 4, synthetic spectra of models with p = 0.2 (model 
A20p0.2) are shown. At t = 20 days (the synthetic spectra are 
similarly scaled, upper panel of Fig. 4), although the absorp- 

7 In Type Ic SNe, R-T instabilities could develop at the (C+0)/(Si+Ni) 
interface (Kifonidis et al. 2000). 
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Fig. 4. — The spectra computed with model A20 and A20p0.2 compared 
with the spectra of SN 1998bw. The synthetic spectra are reddened with 
E(B-V) = 0.032 and scaled assuming fi = 32.76. Upper: Synthetic spectra 
at 20 days after the explosion (shifted upward by 2.0 and 1.0 x 10 -14 , respec- 
tively) and the spectrum of SN 1998bw 18 days after GRB 980425. The flux 
is reduced by a factor of 1.25 after the computation. The vertical lines are 
drawn at 4200, 4700, 5900 and 7500 A. The horizontal lines show the same 
flux level (8.0 x 10 -15 ). Lower. Synthetic spectra at 30 days after the explo- 
sion (shifted upward by 1.6 and 0.8 xlO -14 , respectively) and SN 1998bw 
28 days after GRB 980425. The vertical lines are drawn at 4250, 4800, 6050 
and 7500 A. The horizontal lines show the same flux level (4.0 x 10 -15 ). 



multi-dimensional hydrodynamic and nucleo synthetic mod- 
els, which are computed ab initio, there are not many param- 
eters that we can freely control. The agreement could be im- 
proved by additional modifications, but they would be com- 
plicated and are beyond the scope of this work. Note that our 
input luminosity is higher than that of SN 1998bw at t = 20 
days, making the ejecta temperature too high. If a luminos- 
ity comparable to that of SN 1998bw is used, neglecting the 
possible displacement of the photosphere, the O I-Ca II fea- 
ture becomes as strong as the observation. Even in this case, 
however, mixing with p ~ 0.1 is still required. 



6. CONCLUSIONS 

We have presented the first detailed simulations of the early 
phase optical spectra of realistic jet-like hypernova models. 
The emergent spectrum is different for different viewing an- 
gles. The spectral properties are determined by the combi- 
nation of aspherical abundances and anisotropic ionization 
states. 

Considering the complexity of the problem, the synthetic 
polar- viewed spectra are in reasonable agreement with those 
of SN 1998bw, when ~ 10-20% of the synthesized mate- 
rial is mixed out to higher velocities. The kinetic energy of 
an aspherical model that reproduces, at least qualitatively, the 
spectra of SN 1998bw (E 5i = 20) can be smaller than that of 
a well-fitting spherical model (E$i = 50). This is consistent 
with previous results obtained from models of the LC and the 
late time spectra (Maeda et al. 2006a & c). 

The simulations enable us to predict the properties of the 
early phase spectra of hypernovae viewed off-axis. Compared 
with the spectra seen from the polar direction (red), the spec- 
tra seen from the equatorial direction (blue) have (1) a slightly 
lower absorption velocity, (2) stronger peaks around 4000 and 
4500 A and (3) a stronger Na I A5890 line than in SN 1998bw. 
SNe similar to SN 1998bw in ejecta mass, kinetic energy and 
56 Ni mass should always show spectral features of "hyper- 
novae" or "broad line supernovae", but with some differences 
as described above, depending on the viewing angle. 



tion line at 7000 - 8000 A is still weaker than the observa- 
tion, the Fe lines are strong, causing the broad absorption at 
< 5000 A followed by the peak at - 5200 A. At t = 30 days, 
the strong Ca II feature contributes to the broad absorption at 
7000 - 8000 A, which is comparable to that in SN 1998bw. 

Although the agreement of the spectra is far from per- 
fect, the spectra of the model with mixing are in qualita- 
tive agreement with those of SN 1998bw. Since we use 
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